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Abstract 

A new expression for the Green's function of a finite one-dimensional lattice with nearest neighbor 
interaction is derived via discrete Fourier transform. Solution of the Heisenberg spin chain with 
periodic and open boundary conditions is considered as an example. Comparison to Bethe ansatz 
clarifies the relation between the two approaches. 
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I. INTRODUCTION 



Lattice Green functions have a broad spectrum of applications in physics and allow to 

ly the lattice is assumed to be 



Often the range of relevant 



study dynamic and statistical properties. Usua 
infinite while distances remain finite and discrete, 
excitations allows to consider the system as continuous and apply field theoretic methods. 
8| Accordingly, the integral Fourier transform or infinite Fourier series are commonly used 
as powerful tools for such problems. It therefore appears natural to employ the discrete, or 
finite, Fourier transform (DFT) when the finite size is relevant. However, DFT has mainly 
served for numerical analysis or in fast Fourier transform calculations and much less as 
an analytic method. Most of the known exact results on finite systems have been obtained 



in the Bethe ansatz approach [lo| (see, e.g., 11, 12] for recent review) that is based on a 
different strategy. In 3( the Bethe ansatz solution is compared to the solution by Fourier 
series in the thermodynamic limit of the periodic quantum spin chain. However it is noted 
that the relation between them remains obscure. For instance, the bound two-magnon 
excitation is characterized by the same eigenenergy, but only one eigenstate is found in the 
Fourier approach, instead of the two solutions in Bethe ansatz. This discrepancy has been 



resolved in [13| by an approximate treatment of finite size corrections. In the present paper 
it is shown that the DFT approach recovers the exact solution. The derivation presents 
interest in itself, in view of further extensions. As another example, some results on open 
boundary conditions obtained by Bethe ansatz [ijl are recovered. 

The finite lattice sums are first expressed in terms of infinite Fourier series. For the 
nearest neighbor interaction the series represent an expansion in modified Bessel functions. 
Its Laplace transform results in analytic expressions for the finite lattice Green functions 
and allows to determine the corresponding discrete Fourier transform. 

The standard definitions of the Fourier series for a function periodically continued from 
the interval x G [0; 2n) are 

oo 

/ (z; x) = ^2 c n exp (mi) , (1) 



n=— oo 

2tt 



/ 0; y) exp (-iny) dy, 



1 

Cn ~^Jo 

where the parameters z and 5 (below) take real values. The Fourier expansion of the lattice 
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sum 



m=0 



can then be written as 

oo 

c n exp (m5) 



1 \ -s / n 

— } exp i2nm— 
N \ N 



m=0 



CfcArexp (iSkN) . 



k=—oo 



It has been taken into account that the sum in the square brackets above gives either 1 or 
0, provided n/N is an integer (k) or not. We note that the same property is essential for 
the definition of the DFT of a function a (X) defined on a set of N consecutive integers 
(0,1,.JV-1): 

a W = n S b [ q = AT m ) 6XP {nX) ■ (2) 

m=0 ^ ' 

Then, with the expression of the Fourier coefficients c in (CQ), our lattice sum becomes 

1 N — 1 OO 1 „27T 

- £ / [z- 6 + 2vr^) = E ^ / ex p ( 5 - *)) d y- ( 3 ) 



m=0 fc=— oo 

Now the problem of finite summation has been reduced to integration and series summation. 

To motivate a particular choice of / (z; y) in ([3]) we note that dynamics of different models 
on a finite d-dimensional lattice can be described in terms of Green functions of a generic 
form 

(4) 



G(p;X,...,Y) = 1 Lj2-J2 ^ {lQlXl + '■' + tQdYd) 



N o"o P-e(Q d ,...,Q d ) 
where p is the spectral parameter, e is the lattice dispersion and (X 1: X d ) are the discrete 
coordinates on a periodic lattice (X £ [0, 1, N — 1]). The conjugate " quasimomenta" 
(Qi, Qd) defined as 

2tt A 

= iv mi + iv' mi = "' ' ( 5 ) 

are already present in the l.h.s. of ([3]) and their definition is explained in the example below. 
For simple lattices the dispersion is a sum of cosines 



e = «i cosQi + ... + ad cosQ 



(6) 



Therefore the Green function can be related to a Laplace transform of 



/ (z;y) = exp (z (a x cosyi + ... + a d cosy d )) . 



(7) 
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II. LATTICE GREEN'S FUNCTION 



Let us consider the consequences of d3J) for the choice (JTJ) of / (z; y) corresponding to 
the chain with nearest neighbor interaction. Integration gives the modified Bessel function 
IkN (z) an d ([2]) becomes 

^ iV— 1 00 00 

— ^2 ex p ( z cos Qm) = E IkN ^ exp ( ik ^ = J o w + 2 E IkN ( z ) cos • ( 8 ) 

m=0 k=— 00 k=0 

The parameters z and A can now be continued to complex values, e.g., rewriting Eq. (jSJ) in 



an equivalent form, we obtain a generalization of the Jacobi expansion 

N-l 



15|: 



E ex p (\ (- ex p (t^r) + h exp (-^r) ) ) = E ^ w ^. (9) 

m=0 ^ ^ ^ ' ^ /// — 00 



One can check that the expressions 9.6.33-9.6.40 in 16j are particular cases of (181) or (191) for 



N = 1,2. For instance, taking A = 9 + | in (jSJ) one finds 



~^E exp r COS (^ + 277 + 27r 9) ) = E W*)e*P ex p(^) ( 10 ) 

m=0 ^ ^ ' ' fc=-oo ^ ' 

oo oo 

= I (z) + 2j2 ("I)" hkN (z) cos (2fe0) + 2 ^ (-l) fc / (2fe+1)JV (2) sin ((2k + 1)9). 

k=l k=0 

It is also not difficult to obtain a further generalization of these results by including the 
space dependent factor of the Green function (J3J) and a dispersion of the form: 

e = bcosQ + csinQ, (11) 
so that the propagator can be written as 

m=0 

It corresponds to the Laplace transform of the lattice sum 

L = ^ E ex p (* cos + 27 ^ + ^) ) ex p (< + 2 -i) x ) • ( 13 ) 

By denoting £ = A + r/iV we have 

L = exp (-i V X) 1 £ exp cos (| + 2vr^) ) exp (z (| + 2tt^ x) . (14) 
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Then substitute (JUJ) into ([3]) 

L = exp (— 277X) exp (z/c£) — / exp (z cos y) exp (— z (&X — X) y) dy, 
and finally obtain 

00 

= exp(-^X) IkN+x (z) exp (-ik (A + r]N)) . (16) 

fc=— 00 

The above expressions demonstrate a direct relation of the generalized Jacobi expansion to 
the finite lattice problem. However their main advantage becomes clear by realizing that the 
increasing orders of the Bessel function in these series correspond to the consecutive terms of 
the asymptotic X— expansion of the lattice propagators. I.e., the term with k = describes 
the thermodynamic limit ( X — > 00 ) and the first finite-X correction is contained in the 
term with k — 1. One can also evaluate the convergence of these series, that depends on 
lattice and on the spectral parameter p. We will provide specific examples and explain the 
choice of " quasimomentum" Q and its apparent conflict with periodicity of the lattice. 

The lattice Green function ffl2|) can now be obtained from the Laplace transform of ffTB"]) 
with respect to z 17] 

exp (— 1 77. 1 v) 



/•oo 

/ exp (-pz) I n (z) dz 
Jo 



sinhf 

where n is an integer and we have introduced the notation p = cosh v ( convergence of the 
Laplace transform requires Re p > 1 ). Thus 



sinhv^ 1 exp (z (§ + 2vr^) X) 



X cosh v- cos (# + + 

00 

= S2 exp(-v\kN + X\)exp(-ik(A + r)N) - irjX) , 
fe=— 00 

leads to the following expression for the propagator (fT2l) : 

1_ exp (z (| + 2tt^) X) = exp {vq {N/2 - X M ) + z (M + 1/2) A) 

N ~' coshv ~ cos {§ + 27r f + v) 2sinh'u 

/ exp(t;(X/2-X M )) exp (-y (X/2 - X M )) \ 

Vsinh((wX + i(A + r/X))/2) sinh ((vN - i (A + rjN)) /2)J ' 1 J 
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The integer M in the (fTTI) defines a translation relating an arbitrary X to Xm from the 
main interval: 

X M = 0,1,..JV-1, (18) 
X = X M + MN. 

We note the transformation properties following from the above definition, 

X ^ X + X : X M ^ X M , M — M • 1: 

X -> -X : X M ^iV- X M , M -> - (M + 1) ; 

X->X-X: X M ^ N - X M , M -»• -M. (19) 

After the r.h.s. of (ITT)) has been found, it is now easy to verify that the inverse DFT (i.e. 
taking the summation in X ) indeed reproduces the function in the l.h.s.With ( [191 one also 
finds another useful form of DFT (117)1 : 

2 sinh?; ^ cos((f + 27rf)X) _ cosh ((v + irj) (f - X M ) + i (M + |) A) 
X ^ cosht; - cos (f + 27rf + 77)" sinh (f (v + 177) + if) 

| cosh ((y - ig) (f - X M ) - i (M + |) A) 
sinh (f (f — ^77) — if) 

2zsinh^ sin((f + 27rf)X) _ sinh (fo + irj) (f - X M ) + i (M + |) A) 

X ^ coshw-cos(f + r/ + 2vr^) ~ sinh ({v + i^ f + iA/2) 

_ sinh((^-^) (f -Xm)-^(M+1) A) 

sinh ((« - 177) f - iA/2) ' 1 ' 

The condition cosh v ^ 1 in our derivation was due to the use of the Laplace transform. It 
is however important to mention that the above expressions can be analytically continued to 
arbitrary complex values of the parameters ( v,rj, A ), since the resonance conditions (zeros 
of the denominator in the l.h.s.) correspond to simple poles. For instance, a physically 
important possibility is the region Re p < 1 (p 7^ cosQ), where respective solutions of the 
Schrodinger equation with Im p = are known as scattered states. 



III. THE HEISENBERG RING 



The model of a cyclic chain of X quantum spins S = 1/2 wi p ; 
neighbor interaction was the first to be solved by Bethe ansatz 



h ferromagnetic nearest 
lo| and it offers a clear 



test for the above formal results. The two-magnon subspace of the eigenf unctions of the 
Hamiltonian 

N 

ff=-jX)S < -S i , (21) 

(ij) 

contains the main features of the interaction problem and is described by the amplitude of 
two spins flipped on the sites niand n 2 : 

l<ni<ri2<N 

Within the Bethe ansatz approach a phase shift parameter is introduced to account for 
nearest neighbor interaction in terms of a boundary condition problem. I.e., by matching 
this free parameter in the noninteracting particle form of the ansatz wave function at large 
separation to satisfy the Schrodinger equation when the two overturned spins are nearest 
neighbors. [U] An alternative approach is based on a representation of the Schrodinger 
equation in a form consisting of two contributions: free motion and interaction terms. [3J 

[E - 2J) a{X) + J cos f^j (a{X + 1) + a{X - 1)) = (22) 



J 



PX 

cos ( ) o(0) - a{X) 



{Sx, i + 5 X , n-i) , 



2 

where a(X) represents the amplitude of the relative position X of the flipped spins, E is the 
excitation energy over the ground state with all spins parallel and P is the total momentum 
associated to translation symmetry. 

A(n u n 2 ) = e -^^-a(X); (23) 



R= n 1 ^_n 1 _ X = (n 2 - ni ) = l,...,N-l. (24) 

The split form of fT22|) is convenient to set up a perturbation theory and is widely used in 
condensed matter, when the exact solution is not known. Then the solution is sought in 
terms of Fourier transform. The negative values of X correspond to transposition of the two 
overturned spins, i.e., to the same state. At the same time it is a common practice to chose 
the interval X symmetrically with respect to zero (—N/2, (N — 1) /2 ) and interpreting 
negative values as "left neighbor" and positive values as "right neighbor". The respective 
amplitudes do not necessarily coincide. Therefore here the definition in (Tl8l) is preferred to 
avoid this ambiguity and for a straightforward comparison to Bethe ansatz that follows the 



same counting convention. For instance, if the "right neighbor" is defined by X , then the 
"left" one corresponds to iV — X. 

First we note that quantization of the total momentum P = 2irk/N, k = 0, 1, .., N — 1 in 
( 12"5|) follows from translation of the two spin complex as a whole A(n\, n 2 ) = A(n\ + N,n 2 + 
N). Respectively, the quasimomenta Q in ([T2]) . that are conjugate to relative distance X, 
should be determined by the boundary conditions on the displacement of separate spin flips. 
Thus the amplitude of relative motion in (1231) should satisfy the condition 

a (X) = exp (ink) a{X + N), (25) 

following from the translation A(ni,n 2 ) = A(ni,ri2 + N). According to f[2"5"j) the relative 
amplitude is periodic either on a length of the chain or on a double length, 2N, depending 
on the parity of the total momentum quantum number k = PN/2ir. In the latter case the 
amplitude changes sign after completing the first cycle and may be denoted as antisymmetric, 
a, to distinguish from the former, symmetric mode s, which completes the period after the 
first cycle. Together with the property a (X) = a {—X) due to spin transposition symmetry, 
the above relations define the continuation of the solutions of (1221) from the main interval 
( |24|) to arbitrary values of X. By applying these relations to the Fourier expansion 

Q 

we obtain the condition on Q 

expAiv^±gn =i. 

It determines the two sequences corresponding to the two types of modes 

Qs = lT ] * = 0,1, 

Qa = W + lv ] t = Q > 1 >-> N - 1 - ( 27 ) 

These quasimomenta correspond to rj — and A s = 0, A a = 7r in the general expressions 
([5]). The state counting convention for a one-dimensional lattice is to require that the second 
argument in the amplitude A(n\ ) n2) is larger then the first. This leads to an additional 
relation A(rii,n2) = A(n2,rii + N) imposing a constraint on the amplitude for the main 
interval 

a (X) = exp (ink) a(N-X). (28) 
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For instance, 

a a (X) = -a a (N-X) = -a a (N + X) = a a (-X) . (29) 

One can check that (1281) is automatically satisfied for the amplitude (12T)]) and (1271) . 
Solving ([22]) for the b (Q) in (j2B) . we find 

XQ)=C(P) h C05W ' , 
cosh v — cos Q 

where the constant C (P) is 



cos I — ) — cos Q 



Q' 

The parameter v is related to the eigenenergy E = 2J(1 — cos (P/2) cosh t>) and is deter- 
mined by the compatibility equation 

1 y> COS ((J') [cOsQ' — COS (j)] 

iVcos(^) z V cosh t> — cos Q' 

v 2 y Q 

From the DFT formulas ffT7|) . (TSUI) and (I2T1) we find the eigenenergy equation for each type 
of eigenstate. For instance, 

cos(QJ 1 sinh(o a (Af-2)/2) 



N ^ 



N l—j coshf a — cosQ a 2 sinh (t> a ) cosh {v a N/2) 



After simple algebra we finally obtain 



coshf ajS - cos ( ^ j j cot h ( Nla,S 2 l ^ a ' S ) = sinhf aiS . (30) 



Eqs. (130]) indeed coincide with the Bethe ansatz result (lOj, [Uj. Solutions with real values 
of v , describing the two bound states mentioned earlier, merge in the thermodynamic limit 
to: v = — In (cosP/2), E = J sin 2 (P/2) . At any given P equation (1301) has a number ~ 
of solutions with imaginary values of v, or Re p < 1, which correspond to scattered states. 
These solutions are characterized by an oscillating space dependence, as is clearly seen from 
( 120]) (see also (|3T]) below), unlike the smoothly decaying dependence of the bound states. At 
finite iV they form two distinct classes of energy bands, shifted with respect to each other 
by ~ J/N. All the solutions can be cast in a unified form by using the relation of A to the 
total momentum (A = PN/2 ): 

a (X) ~ cosh (v - \X\\ + i^j , (31) 



where v = v (P) corresponds to a particular solution of (1301) at fixed total momentum P. 
One can see that the phase shift of tt/2 ( stemming from PN/A ) between nearest values 
of P survives in the thermodynamic limit. This distinction in symmetry of eigenstates was 
missed in the previous Fourier series analysis 0, 18] because the effect is contained in the 



phase shifts of quasimomenta Q (1271) that vanishes in the thermodynamic limit. 

IV. OPEN BOUNDARY CONDITIONS 

To consider another illustration of the formalism derived above we will include, according 
to 14|, also arbitrary end fields (/i and v) into the anisotropic Hamiltonian on a chain of 
length L ( i,j = 0, ...,L — 1 ). Only a single magnon excitation will be considered, as the 
interaction is already present due to noncyclic boundary conditions: 

H = ~ E (W^ + i (S+S7 +1 + S7S+ +1 ) + ^ + v si_^j ■ 

The parameters are dimensionless, i.e., scaled with the exchange constant J x , and the energy 
of the state with all spins parallel is E = —r) (L — 1) /2 — + v) /2. The Schrodinger 
equation is represented in the form similar to ( 1221) with the "interaction" on the r.h.s.: 

{E-E )a{j) 



Va (j) - (j - 1) - -a (j + 1) 



= S j=L ^ Qa (L) - a (L - 1) Q77 - 1/ 

+ *i=oQa(-l)-o(0)^-A*)). (32) 

As can be seen, the nonphysical amplitudes a (—1) and a (L) cancel each other and in the 
DFT approach are in fact defined by the Eq. (132]) and the transformation properties (Tl9l) . 
For instance, these amplitudes may differ from Bethe ansatz solution, as is the case of the 
Heisenberg ring. The two approaches, however, give the same solution for the physical 
interval. In DFT we need to consider an interval in X which can be periodically continued. 
However unlike the cyclic chain, here the last spin in the open chain is not at the same time 
the "left" nearest neighbor of the first spin, i.e., there can be no periodicity transformation: 
X — > X + L. Instead, our chain can be viewed as a half of the cyclic chain (ring) of length 
N = 2L ( j = 0, .., 2L - 1) where the two links (2L - 1, 0) and (L - 1, L) are broken. With 
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this definition we can apply the DFT formulas (j2J) and (jTTl) to solve (|32|) . IN this case 
A = rj = and the quasimomentum is obviously defined as Q = 27rm/N, m — 0, 2L — 1. 
An immediate consequence of transformation properties (I19p is that a (— 1) = a (1) , while 
a (L) has to be determined from (1^21) . The Fourier amplitude resulting from Schrodinger 
equation is 

(cosh v — cos Q) (cosh u — cos Q) 

where coshw = r] — (E — E ) and the two constants are 

d = \a (-1) - a (0) Q^-/i 

C 2 =~a(L)-a(L-l) Qij - ^ . 

Substituting the amplitude 6 (Q) into the definition of the constants one obtains a linear 
homogeneous system of equations 

N-l 



[exp (-iq) - (rj- 2//)] 



-c 2 — y 



2N ^— ' cosh f — cos (/ 

m=0 ^ 

AT-1 



exp (-ig (L - 1)) [exp (-zg) - (77 - 2/i)] 



' 2N ' cosh t; — cos o 

m=0 ^ 



C Jl + _Ly [ ex P^)-^- 2 ^ ) (34) 



2iV ^— ' cosh f — cos o 

m=0 ' 



= _ C r J_ y^ 1 exp (zg (L - 1)) [exp (ig) - (77 - 2u)} 
2N ^— ' cosh v — cos q 

m=0 ^ 

which defines the eigenenergies and eigenf unctions of the equation (|32|) . From (fT7|) we have 



N-l 



1 E 



exp (— ig) cosh (f (L — 1)) 



iV ^— ' cosh f — cos sinh t> sinh (t> L) 



m=0 
jV-1 



cosh (vL) 



iV ' cosh v — cos sinh v sinh (t> L) 

m=0 ^ v y 

1 ^-4 exp (— iq (L — 1)) cosh (t>) 
AT 2^ 



iV ' cosh u — cos sinh w sinh (t>L) 

m=0 ^ v ' 

1 exp (— iq (L — 2)) cosh (2t>) 

N ' cosh v — cos sinh t> sinh (t> L) 

m=0 ^ V 1 
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With these expressions we obtain from (I34|) the eigenenergy equation of 



e 2v { L-x) = (e- -( V - 2^)) (e- - (77 - 21/)) 

(e w - (77 - 2/x)) (e 1 ' - (77 - 2v)) ' 1 J 

The equation is invariant under v «-> — 7j or /i «-» z/. By taking the ratio C1/C2 from (I34|) one 
obtains from ffTTl) the space dependence of the amplitude: 

a (X) = 7 [sinh (v (x + 1)) — (77 — 2/j) sinh (t>x)] . (36) 

The factor 7 = 7 (u, L) can be found from normalization. It is now easy to checked that 
(|36l) is equivalent to Bethe ansatz result provided t> satisfies the eigenenergy equation (|35l) . 
Interchanging the end fields /i <-> v induces an obvious transformation (up to a phase factor) 
a (X) — > a (L — 1 — X) which amounts to changing the direction of counting the sites. As 
noted in, 14( the solution is particularly simple for the condition 

(77-2/i) (T7-2Z/) = 1, 

when one gets e~ 2v in the r.h.s. of (|35|) and consequently: 

71m 

v = i— — , m — 0, 2L — 1. 

The wave function represents a superposition of two harmonic waves. Looking at the energy 
of the excitation 

E — E = 77 — cosp, 

one can immediately see that in the case of anisotropic ferromagnetic interaction with 77 < 1 
(J z < J x ) the ground state is unstable against long wavelength magnons. The harmonic 
(purely oscillating) solution is also obtained in the limit of strong end fields fi, v — > 00 
(equivalent to fixed boundaries) since (135]) becomes 

e 2v(L-i) = 1 . v = ipi p= m = o, 2L - 3. 

Lj — 1 

Due to the symmetry of the chain mentioned above, only half of solutions is linearly inde- 
pendent. An additional solution corresponds to bound state recovering the total number of 
modes L. For an arbitrary choice of parameters the solution will consist of scattered waves 
( Ret> = 0) and bound states ( Imv = ). 
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V. CONCLUSIONS 



A new expression for the finite ID lattice Green function has been derived by the discrete 
Fourier transform approach, by establishing a generalization of the Jacobi expansion. We 
note its qualitative advantage compared to the approach based in infinite Fourier series that 
it make accessible a new region of parameters, not available in the thermodynamic limit, 
but having an important physical meaning (e.g., scattering wave solutions). In the Bethe 
ansatz approach the phase 9 is introduced into the assumed form of the wave function and is 
then determined by substitution into the Schrodinger equation. As a result, the Bethe phase 
becomes a function of momentum (and consequently of energy) and contains all the informa- 
tion on the interaction in the system. In DFT the sequence is just opposite: no assumption 
is made on the wave function. The eigenenergy (Schrodinger) equation is first solved in 
the momentum space and then the wave function is obtained via Fourier transform into the 
direct lattice space. Then the identification of the Bethe phase is achieved by comparing the 
two functions, 9 = ivN + 7tm. The DFT allows to describe the scattering states on the same 
footing as the bound states, and at the same time to keep a physically transparent structure 
discussed above. It is known that the bound states are well separated from continuum of 
scattered states in the thermodynamic limit of (l2Tj) . 181 ] while the description of scattered 
states requires special treatment in a Lipmann-Schwinger type approach jj. In the DFT 
approach it becomes clear that the qualitative change of the wave function is controlled 
by the effective magnon-magnon interaction described by the parameter v (e.g., the r.h.s. 
of Eq. (1221) ): the bound states are formed when this interaction is attractive, while the 
scattered states are due to repulsive effective interaction (see ll3|). One can show that the 
critical line separating the two classes of states corresponds to the vanishing of interaction 
(i.e. the magnons are free if their dispersion crosses this line) and coincides with the lower 
boundary of the continuum of scattered states. The wave function of the antisymmetric 
solutions has nodes on the direct lattice and therefore the corresponding states are more 
loosely bound. At sufficiently long wavelength of the excitation the bound states are very 
close to scattered states and the antisymmetric ones can indeed become unstable, i.e., decay 
into scattered states. This mechanism is responsible for the so-called non-string behavior of 
the solutions of the Bethe equations at finite N. [19( Thus also the physical interpretation 
of the Bethe solutions becomes more transparent. The same approach allows to reproduce 
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the BA solution for an open chain, which is amenable to Fourier expansion by doubling the 
length of the chain. This can be viewed as a cyclic ring with broken links (or impurities) 
introducing a scattering of magnon excitations and decay of the wave function with distance. 
The approach described in the paper presents interest for further development, since by its 
construction it is not limited to one space dimension. 
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